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The origin of pseudobulges in cosmological simulations of 
galaxy formation 
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ABSTRACT 

More than half of nearby disc galaxies have pseudobulges, instead of classical bulges 
that are though to be end-products of galaxy mergers. Pseudobulges are presumed to 
develop over time as a result of secular evolution of galaxy discs. We report simulations 
of galaxy formation, in which two disc galaxies with disky pseudobulges have formed. 
Based on the profile decomposition, the bulge-to-total mass ratio of the simulated 
galaxies is 0.6 for one galaxy and 0.3 for the other. We find that the main formation 
mechanisin of the pseudobulges in our simulations is not the secular evolution of 
discs but high-redshift starbursts. The progenitors of the pseudobulges form as high- 
redshift discs with small scale lengths by rapid supply of low angular momentum 
gas. By redshift 2, before the main disc formation, pseudobulge formation has largely 
completed in terms of mass. The secular evolution such as bar instability accounts for 
about 30 % of the bulge mass for one galaxy and only ^ 13 % for the other but does 
affect the final shape and kinematic properties of the pseudobulges. 

Key words: methods: numerical - galaxies: formation - galaxies: bulges. 



1 INTRODUCTION 

Pseudobulges are characterised by distinctive features 
i-H ■ from class ical bulges that are though t to be merger 
L^ ■ remn a nts (iKormendy fc KennicuttI |2004| : iNaab fc Truiilld 
, ^ ! I2OO6I : iHopkins et al.l I2OI0I '). Pseudobulges often show 
substantial rotational support and disky or boxy/peanut 
edge-on shapes. These features support the idea that 
'Seudobulges have formed vi a secular evolution of discs 



t 



, pseudobulges have tormed vi a secular evolution 01 discs 
(jKormendv fc KennicuttI 120041 ') caused by non-asymmetric 
structure in a galaxy disc s uch a s spiral arms and a 



form via hierarchical merging that naturally produces 
classical bulges. Observationally however more than 
half of b ulges of nearby la r ge di s c galaxies are pseu- 
dobulges (iGraham fc WorlevI l2008l : IWeinzirl et all l2009l : 
iKormendv et al.l uOlOT ). This observation al fact could be 
a ch a llenge to the standard cosmology (|Kormendv et al.l 
I2OIOI : [Peebles fc Nusseij [2OI0I ) . Contrarily to the intuition, 
cosmological simulations of Milky Way-sized galaxy forma- 
tion often find that bulges show surface bri ghtness profiles 
with t he Sersic indices small er than 2 (lOkamoto et al.l 



bar ( e.g. Combes fc Sanderg 1981 



199d: ICombes fc Ehnegreenl Il993l : 



Pfenniger fc NormarJ 



iDebattista et al . ^ 
Athanassoulal l2005l l o r by accretion of clumps forme d 



2004 



20051: iGovernato et al] l2007l : IScannapieco et all 12011 . 
Guedes et al.l I2OIII): such bulges ar e usua ll y classified as 
pseud obulges (JFisher fc Dror4 12008|, I2OIOI : iweinzirl et"al] 



gas rich unstable disc (JNoguchill 19991 : llnoue fc Saitohll2012h . 
We call evolution that is originated by internal processes 
'secular evolution' in this paper regardless of its timescale in 
order to distinguish evolution driven by external processes 
such as mergerqj. 

In the standard picture of galaxy formation, galaxies 
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^ Indeed, secular e volution can have timescale mu ch shorter than 
a Hubble time ( see iKormendv fc KennicuttI |2004| . for a review). 
Ilnoue fc Saitohl 1I2OI2I ) showed that the formation timescale of a 
clump origin pseudobulge is rather short and they distinguished 
it from secular evolution. 



[2OO9I). How these bulges form over the cosmic time is 
however remain unclear. 

Investigating formation processes of the simulated 
bulges of disc galaxies should provide important clues 
to understand why pseudobulges are so common both 
in the local Universe and simulated universes. Simula- 
tions for this purpose must resolve detailed structure of 
galaxies, such as shape of bulges. Only recently cosmo- 
logical simulations with sufficient r esolution have become 
possible LO kamot o fc FrenkI l2009l : lOkamoto era!] I2OIOI : 
iGuedes et al.ll201ll ). 

The simulations presented in this paper include a 
number of baryonic processes known to be relevant to 
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galaxy formation, such as radiative gas cooling, star for- 
mation, sup ernova feedback , and chem ical evolution. As 
descri bed in lOkamoto et al.1 (|2005|) and IScannapieco et alj 
l|201H ). there are many uncertainties in modelling of these 
'sub-grid' physics. In particular, feedback is one of the 
most poorly understood processes, while i t most strongly 
affects the properties of simula ted galaxies (jOkamoto et al.l 
I2OO5I . I2OIOI : ISchave et aLllJOlol ). We make use of a model 
that has already had success in reproducing some prop- 
erties of the Local Group satellite galaxies, including 
the l uminosity funct i on an d luminosity-metallicity rela- 
tion lOkamoto et al.l [2010|) and observed properti es of 
high-redshift Lv man-g emitters (|Shimizu et al.l I2OIII ) and 
sub-mm galaxies (|Shimizu et al.ll2012r ). 

The paper is organised as follows. In Section 2, we de- 
scribe our simulations, providing brief descriptions of our 
modelling of baryonic processes. In Section 3, we first present 
properties of simulated bulges at redshift to show that they 
have the characteristics of pseudobulges. We then show red- 
shift evolution of the simulated galaxies and their bulges in 
order to address formation process of the pseudobulges. We 
compare our results with observations in Section 4. Finally, 
in Section 5, we discuss the results and summarise our main 
conclusions. 



e = 1.03 kpc until redshift 3, and then fixed in physical coor- 
dinates as e = 0.257 kpc in both simulations. We adopt the 
same softening length for the high-resolution dark matter, 
gas, and star particles. 



2.1 The simulation code 

The simulation code is based on an early version of a 
Tree-PM SPH code, g adget-3 and it is e xactly the same 
as the one used in IParrv et al.l ((201^). The baryonic 
processes are modelled so as to reprod uce the properties 
of th e Local Group satellite galaxies (|Okamoto fc FrenkI 
I2OO9I : lOkamoto et al.1 I2OI0I I. The simulations include 
metallicity dependent radiative cooling and photoheating 
(jWiersma et al.l 12003 ) in the presence o f a time-evolving, 
spati ally unif orm ultraviolet backgr ound (jHaardt fc Madaru 
l200ll . see also lOkamoto et al.ll2008al ). The star formation oc- 
curs in dense gas with wh > 7^th, where nth is the threshold 
gas density above which the star formation is enabled. The 
star form ation efficiency i s set to reproduc e the Kennicutt 
relation (jKennicuttl 119981 ) as des cribed in Okamoto et al.l 
(|2008bl ) (see also G3-T0 model in lScannapieco et al.ll2012l '). 
The simulations include both Type II and la supernovae 
and follow the evolution of the chemical elements produced 
by these two types of supernovae and by AGB stars. 



2 THE SIMULATIONS 

To investigate the properties of Milky Way-sized 

galaxies, we select two of the six haloes from the 

Aquarius project (|Spring:el et al.l I2OO8I ): 'Aq-C' and 
'Aq-D' in their labelling system. We hence assume 
a ACDM cosmology with the following parameters: 
f^o = 0.25, ^A = 0.75, Qb = 0.045, erg = 0.9, n^ = 0.9, 
and a Hubble constant of Ho — 100 h km s~^ Mpc~^, 
where h = 0.73. These parameters are consis tent with the 
WMAP 1- and 5-year result s at the 3 a level (jSpergel et al.l 
I2OO3I : iKomatsu et al.l [2009|). These haloes are extracted 
from a cosmological periodic box of the side length of 
100 h^^ Mpc, and were chosen to have masses close to that 
of the Milky Way's dark matter halo (~ 10^^ Mq). These 
haloes are relatively isolated and do not h ave neighbours 
excee ding half its mass within 1 h~^ Mpc (JNavarro et al.l 
l2010t ). 

As in Aquarius, we use the resimulation technique; the 
initial density field of the parent simulation is recreated 
by adding appropriate additional short wavelength per- 
turbations in the region out of which the halo of interest 
forms. In this region, we also place gas particles which are 
used to perform Smoothed Particle Hydrodynamic (SPH) 
calculation. The region external to this is populated with 
high-mass dark matter particles, the function of which is to 
reproduce the appropriate tidal fields. 

The high-resolution dark matter particle mass is 2.6 x 
lO'^ Mq and 2.2 x 10^ Mq, respectively for Aq-C and Aq-D, 
and the original SPH particle mass is thus 5.8 x lO** Mq 
and 4.8 x 10* Mq. Note that the SPH particle mass can 
change owing to star formation and feedback. Gravitational 
softening length is kept constant in comoving coordinates as 
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feedback 

winds 

Pavel 



f rom supernovae ex p losion 




particle 



receives an amount of energy AE from Type II supernovae 
during a timestep. At, this particle is added to winds with 
a probability, pw ~ AE / {^mga,Bvl,) , where TTigas is the 
mass of the SPH particle and Ww is the initial wind speed. 
This is given as v„ — 5a, where a is the one-dimensional 
velocity dispersion of the dark m atter particles aroun d 
the gas particle (vw5cr model in lOkamoto et al.l [2010|). 
The wind particles are decoupled from the hydrodynamic 
calculation for a short time in order to allow them to escape 
the high-density star-forming regions. When the density 
has fallen to nu = 0.01 cm^^, the particles feel the usual 
hydrodynamic force again. If they do not reach sufficiently 
low density after a time 10 kpc/«w, they are recoupled to 
the hydrodynamic interactions anyway. 

This expression of the wind mass-loading implies 
that the wind mass generated by a Type II supernova is 
proportional to a~^, i.e. per unit of star formation, less 
massive galaxies blow more winds than their more massive 
counterparts. This is a key feature to reproduce properties 
of the Local Group satellites, including the shape and 
approximate normalisation of their luminosity function a nd 
the luminosity-metallicity relation (jOkamoto et al.ll2010l ). 

The model for supernovae driven winds u sed h ere 
differs from the original one by lOkamoto et al.l ((20101) in 
two ways. Firstly, we allow all gas particles, not just those 
above the star formation threshold density, to be added 
to the wind if they receive feedback energy. The original 
prescription can result in a variable wind mass-loading 
depending on how well the star-forming region is resolved. 
Secondly, only type II supernovae contribute to the winds. 
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Type la supernova energy is added to the gas as thermal 
energy. 

There are several more changes to the code from the 
one used in the previous studies. Firstly, the star formation 
threshold density is now a function of the numerical resolu- 
tion while it was fixed to nn =0.1 cc~^ in the previous stud- 
ies. For the resolution of the simulations presented in this pa- 
per, we employ nu = 1.6 cc~^ as the star formation thresh- 
old density (see Appendix. \C\ for more details). Secondly, 
we h ave implemented a timestep limiter (jSaitoh fc Making 
12003 ) that reduces the timestep of a gas particle if it is 
too long compared to the neighbouring gas particles. With- 
out this limiter, halo gas particles which have much longer 
timesteps than high-speed wind particles often do not feel 
the force from the wind particles (i.e. do not notice the wind 
particles), while the wind particles feel the drag from the sur- 
rounding halo gas particles. Lastly, we have added artificial 
conducti vity in orde r to capture the instabilities at contact 
surfaces (|Pricdl2p08l) together with the tim e dependent ar- 



tificial viscosity ( Morris fc MonaghanI 



19971): otherwise SPH 



cannot deal with such instabilities ( Okamoto et al.l 120031 : 
lAgertz et al.ll2007l ). All these changes have significantly im- 
proved the numerical conver gence as we show in Appendix 
El (see also iParrv et al.ll2012l ) . 



3 RESULTS 

In this section, we first analyse the properties of the simu- 
lated bulges in order to show that the simulated disc galax- 
ies have pseudobulges. We then show the evolution of the 
galaxies and their bulges to understand pseudobulge for- 
mation. The global properties of the simulated galaxies at 
redshift are summarised in Appendix [X] In the following 
analyses, the ^-direction is chosen to be parallel to the an- 
gular momentum vector of stars within 5 % of the virial 
radiu^J at given redshift. Note that the direction of the an- 
gular momentum vector significantly changes with redshift 
(I Okamoto et al.l l2005l : IScannapieco et al.l l2009l : I Sales et al.l 
[2013) as we will show later. The x-direction is to be parallel 
to the major axis of the surface stellar density distribution 
in face-on. We also define galaxy radius as 10 % of the virial 
radius. 



3.1 Shapes and kinematic properties of the bulges 

For classifying bulges, the Sersic profile fitting for a surface 
stellar density profile is frequently used: 



E(r) — Ee exp 



-b„ 



Re 



1.0 



(1) 



where Re is the effective radius and Ec is the surface density 
at Re, respectively, and n is the Sersic index. The parameter 
b„ is related to n by b„ ~ 2n — 0.32 4. Bulges with n < 2 are 
usually classified as pseudobulges (jKormendv fc KennicuttI 
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Figure 1. Azimuthally averaged surface stellar density profiles. 
The left and right panels show profiles of Aq-C and Aq-D galaxies, 
respectively. The profiles of each galaxy is fitted by a combination 
of the Sersic and exponential profiles. The surface density profiles 
are indicated by the triangles and the best fit profiles are shown 
by the solid lines. The contribution from Sersic bulges (dotted 
lines) and the exponential discs (dashed lines) are also plotted. 
The bulges are well fitted by the Sersic indices, n ~ 1.2 and 1.4, 
for Aq-C and Aq-D, respectively. 



I2OO4I : iFisher fc Dronll2008l . I2OI0I : IWeinzirl et alll2009l fl. In 
Fig. [T] we show the surface density profiles of the simulated 
galaxies from face-on {x-y projection) at redshift 0. We 
fit each profile by a combination of the Sersic bulge and 
exponential disc profiles. We obtain n ~ 1.2 and 1.4 for the 
bulges of Aq-C and Aq-D, respectively. These bulges hence 
have pseudobulge-like profiles. According to the profile 
fitting, the bulge-to-total mass ratios, B/T, are 0.6 for 
Aq-C and 0.3 for Aq-D. These values suggest that Aq-C is 
an early-type disc galaxy and Aq-D is presumably classified 
as an Sb galaxy. The bulge stars dominate mass within 
3 kpc from the galaxy centres; more than 92 % and 86 
% for Aq-C and Aq-D, respectively. We thus regard this 
region as a bulge. 

In FigO we show the face-on (upper panels) and 
edge-on (lower panels) surface stellar density maps. From 
the face-on density maps, it is evident that Aq-C has a 
bar, while Aq-D shows almost axisymmetric surface stellar 
density distribution except in the central region (r < 2 kpc) 
where we see weak non-axisymmetricity. In the edge-on 
density maps, both bulges appear to be quite 'fiat'. This 
is not what one would expect for a classica l bulges; both 
galax ies seem to posses 'disc-lik e' bulges (jAthanassoulal 
|2005| . or 'disky pseudobulges' in lErwinI l2008l 's terminol- 
ogy). A further evidence for a dominant central disk- like 
structure is t he diamond shap ed isodensity contours in 
Aq-D's bulge (|Pohlen et al.ll2000l ). The inner part of Aq-C's 
bulge also shows disky, diamond isodensity contours, while 
the outer part has 'boxy' features (2 < |a::| < 4.5 kpc and 
1 S; kl ^ 2 kpc). The boxy /pe anut shape is related to 
bar structure (jAthanassoulal |2005| ) and therefore the boxy 
isodensity contours are consistent with the fact that Aq-C 



^ The virial radius is c alculated based on the spherical collapse 
model llEke et al.lll996l) 



^ Not e that this classification at n = 2 was questioned bv lGrahaml 
1120111 ) because bulges ranging from n = 0.5 to 10 obey a single 
curved relation in the Re — n plane. 
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Figure 2. Surface stellar density maps of the bulges. The upper 
panels show the face-on views {x-y projection) and the lower pan- 
els show the edge-on views {x-z projection). The contour levels 
are chosen to highlight the bulge shapes, which arc clearly 'disky' 
in the edge-on views. 



has a bar. 

Next we investigate the kinematic properties of the 
simulated bulges. In the top panels of Fig. [S] we show 
the mean line-of-sight velocity maps from edge-on view. 
Both bulges have significant rotation but the shapes of 
the isovelocity contours are rather different between them; 
Aq-C's contours are more vertical than Aq-D's, implying 
the rotation properties of Aq-C's bulge is more cylindrical. 
This is presumably due to the presence of the bar. 

We plot the line-of-sight velocity profiles (rotation 
curves) at the disc plane {\z\ < 0.1 kpc) in the middle panels 
and the line-of-sight velocity dispersion profiles in the bot- 
tom panels. We find that the inner part of Aq-C's bulge is 
rather kinematically hot; the velocity dispersion profiles are 
sharply peaked at a; = 0. The region is supported by the 
velocity dispersion of the stars rather than the rotation. We 
also plot the vertical (z-direction) velocity dispersion of the 
stars in the disc plane. The vertical velocity dispersion is 
only half the line-of-sight one. Therefore the orbits of stars 
in the inner part of the bulge are largely non-circular but 
confined in a thin oblate. On the other hand, Aq-D's bulge is 
kinematically cold. The line-of-sight velocity dispersion pro- 
file is very fia t and it seems typical of the disky pseudobulges 
(JErwinl [2003) • The cold, highly rotationally supported fea- 
ture of the bulge implies that the pseudobulge is a disc-like 
bulge. 



3.2 Evolution of the galaxies and their bulges 

In this subsection we investigate how the galaxies and 
their bulges form and evolve in the simulations. Fig. [4] 
shows stellar and gas distribution in the main progenitors 
of Aq-C and Aq-D galaxies around redshift 3, 2, 1, 0.5, 
and 0. Again we chose the z-direction to be parallel to the 
angular momentum vector and s-direction to the major 




2 ,0,2 

X [kpc] 



Figure 3. Kinematic properties of the bulges from edge- 
on view [x-z projection). Top panels: The mean line-of-sight 
velocity maps. The contour levels correspond to Ktiioa)! = 
0, 25, 50, 75, 100, 150 km s'^. Middle panels: The line-of-site 
velocity profiles at the disc plane (|2;| < 0.1 kpc). Bottom panels: 
The line-of-sight velocity dispersion profiles at the disc plane. We 
also show the vertical (2-direction) velocity dispersion profiles at 
the disc plane by the red dashed lines. 



axis. At redshift 3, both galaxies show irregular morphology, 
characterised by an inner disc and misaligned outer disc 
or ring. The misalignment is caused by the change of the 
direction of angular momentum vector of newly accreted 
gas as we will show later. The components that eventually 
evolve into the discs at redshift start forming around 
redshift 2. Aq-C forms a bar around redshift 1 and the 
bar survives until redshift 0. In Aq-D, a weak (and short) 
bar-like feature can be seen at redshift 2, which disappears 
by redshift 1. Aq-D's late-time evolution below redshift 1 is 
characterised by clumpy star formation. 

Our main result is depicted in Fig. [5] where we show 
the surface stellar density profiles of the main progenitor 
galaxies from face-on view at several redshifts. Redshift 
evolution of the surface density profiles indicates that a 
high-redshift disc with a small scale length forms first, 
then a component that eventually evolves into the disc at 
redshift develops around it. In fact, the Sersic bulges are 
always disc-like {n < 2). It is evident that there has been 
only a little evolution in the bulge mass since redshift 2; 
the bulge mass at redshift 2 already accounts for ~ 70 % of 
that at redshift for Aq-C and ~ 87 % for Aq-D. Because 
the discs seen at redshift start forming around redshift 
2, we conclude that the main channel of the pseudobulge 
formation in our simulations is not secular evolution of the 
main (or outer) discs, although secular evolution seems to 
have non-negligible contribution (~ 30 %) to the mass of 
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Figure 4. Evolution of the galaxies. The side length of each panel is set to 20 % of the virial radius of the host dark matter halo at 
a given redshift. The galaxy radius which is defined as 10 % of the virial radius is indicated by a horizontal line in each panel. The 
left and right concatenated panels, respectively, show Aq-C and Aq-D galaxies. From left to right, the surface stellar density in face-on 
{x-y projection) and in edge-on {x-z), and surface gas density in face-on and in edge-on are presented. Brighter colour is used for higher 
density. From top to bottom, we show the main progenitor at redshift 2.9, 2.0, 0.99, 0.5, and 0. The 2-direction is chosen to be parallel 
to the angular momentum of stars within 5 % of the virial radius at each redshift. The a;-direction is to be parallel to the major axis as 
in Fig. [2] and [3] 
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stars: the stars lie within 3 kpc from the galaxy centres at 
redshift 0. In order to distinguish the stars that are born in 
the main progenitor (in situ) from the stars that are brought 
by accreted sat ellites (accreted), we identify subhaloes by 
using SUBFIND (jSpringel et al.ll2001f ) and construct merger 
trees. We define subhaloes as systems that consists of at 
least 32 self-bound particles. To identify clumps formed by 
the self-gravitational instability of gas (see Aq-D in Fig. Q, 
we also utilise subfind; a group of self-bound particles that 
contains more than 32 particles found within the galaxy 
radius, O.lrvir, whose baryon fraction is more than 90 % is 
regarded as a clump. 



Figure 5. Evolution of the surface stellar density profiles. The 
surface stellar density profiles of Aq-C (left) and Aq-D (right) 
galaxies and their main progenitors obtained from face-on views 
are presented. The profiles around redshift 0, 1, 2, and 3 are indi- 
cated by the black plus signs, blue asterisks, green diamonds, and 
red triangles, respectively. Each profile is fitted by a combination 
of a Sersic bulge and an exponential disk, which is shown by the 
line. The Sersic indices are shown in each panel. 



Aq-C's bulge, which is the one that has the bar at low 
redshift. 

In Fig. [6l we present the formation histories of bulge 



As shown in the top panels of Fig. [6l most of the bulge 
stars are born in starbursts between redshifts 2 and 6 in 
Aq-C and between 2 and 4 in Aq-D. Moreover, most of the 
bulge stars are in situ. This fact implies that at least 'dry 
mergers' do not contribute to the formation of the bulges. 
There is little star formation activity below redshift 1 in 
Aq-C's bulge, while Aq-D's bulge show slight activity at 
low redshift. We find that about half the star formation in 
Aq-D's bulge below redshift 1 occurs in the gas clumps. 
The clumps sink to the gal axy centre by dynami cal friction 
and add mass to the bulge (jlnoue fc Saitohll2012r ). However, 
the contribution from the clumps is only about 10 % and 
not as large as in the 'clump origin pseudobulge' reported 
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Figure 7. Spin flips: Evolution of the orientation of the galax- 
ies. The cosine of the angle between the direction of the spin of 
the galaxy at redshift and that of the main progenitor at each 
redshift is plotted as a function of redshift in each panel. The 
direction is computed using the stellar component within O.Srvir, 
where stars in the satellites are removed from the analysis. 



Figure 6. Distribution of formation times of bulge stars and 
merging histories of galaxies expressed in terms of redshift. Stars 
lie within 3 kpc from the galaxy centre at redshift are identified 
as the bulge stars here. The left and right panels show Aq-C and 
Aq-D galaxies, respectively. The upper and lower panels show the 
same data in differential and cumulative form, respectively. The 
black solid lines indicate all the bulge stars, while the red dashed 
and blue dotted lines respectively indicate those born in the main 
progenitors {in situ) and those born in the satellites (accreted). 
The cyan dot-dashed lines show in situ stars which have formed 
in gas clumps. In the upper panels, we show the merging mass 
ratios by the vertical magenta lines. We show all the mergers 
whose mass ratios are greater than 0.01. 



bv llnoue fc Saitohl (|2012h . 

In order to investigate what triggers the starbursts, 
we study the merging histories of the galaxies. We define 
galaxy mass, Afgai, as the sum of the stellar mass, M«, and 
the interstellar medium (ISM) mass, Mism, in a subhalo, 
were we regard the gas with high density (nu > 0.1 cc~^) 
and low temperature (T < 3 x 10* K) as the ISM. The 
star-forming gas (nn > 1.6 cc~^) is also counted as the 
ISM. The merging mass ratio is defined by the second 
massive progenitor to the main progenitor mass in a merger 
(the vertical lines in the upper panels of Fig. (Gjl . 

We find that the star formation activity is not strongly 
correlated with mergers. We thus conclude that the star- 
bursts are triggered by the rapid supply of gas at high 
redshift instead of the mergers (see Appendix |B]|. We 
also find that both galaxies have not undergone mergers 
with mass ratio greater than 0.1 since redshift 4. We have 
confirmed that the major mergers at high redshift are wet 
(App endix [B]| ; wet merg ers are not as destructive as dry 
ones (|Hopkins et al.luOlOl ') and someti mes form discs by de- 
posit i ng orbital an gular mom entum ( jSpringcl & H ernquistI 
l2005l : [Robertson et al. 200j; iGoverna to ct al. 200i). The 
gas richness at high redshift seems to prevent the galaxies 
from forming classical bulges. 

Finally, we study evolution of the orientation of the 



main progenitor galaxies. During the growth of a dark 
halo, the direction of its spin undergoes rapid changes 
(jBett fc Frenkll2012l ). The direction of the angular momen- 
tum vector of newly accreting gas also chang es dramatically. 
This induces the 'flip' of a disc (lOkamoto c t al. 20 051) and 
some ti mes destroye d a pr eexisting disc (jScannapieco et al.l 
I2OIH ). ISales et al. (120121) po inted out, by using gimic 
simulations ( Grain et al.ll2009l ). that many bulges in Milky 
Way-sized haloes form owing to this misalignment of 
newly-accreted gas. 

In FiglXl we show the evolution of the galaxy orien- 
tation. We define the spin of a disc by using the direction 
of the angular momentum vector of the stellar component 
within 0.5rvir, where we ignore the stars in satellites to 
avoid misdefining the galaxy orientation due to the orbital 
angular momentum of the satellites. The orientation of the 
galaxies rapidly changes at high redshift by the change of 
the direction of the angular momentum vector of newly- 
accreted materials. 

The fiips displayed in Fig. [7] cause to form the mis- 
aligned discs for example seen at redshift 3 in Fig. |3] The 
rapid change of the spin direction at high redshift (> 2) 
prevents disc formation and forms pseudobulges. At redshift 
below 2, at which the host haloes are well established, 
the fiip is subsided and the discs develop around the 
pseudobulges. 



4 COMPARISON WITH OBSERVATIONS 

In this section, we compare the simulated galaxies and 
bulges with observations. We firstly investigate whether 
the simulated galaxies are too bulge-dominated to host 
pseudobulges or not. Galaxies with pseudobulges are 
in general have smaller B/T than those with classical 
bulges (e.g. [Graham fc Worlevll2008l : JFisher fc DrorvlbOOSi : 



Fisher et al.l 120091: IWeinzirl et aJT 120091 '). For example 
Fisher fc Drorvl ( 20081 ) showed that most of the galax- 
ies with pseudobulges have B/T < 0.4 in l/-band and 
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Table 1. Size, mass and star formation rate of the simulated 
bulges at redshift 0. 



Galaxy iJxs 
(kpc) 



Aq-C 
Aq-D 



A^xs 
(Mq) 



Mxs 



Mxs/Mxs 



3.2 

2.7 



2.2 X IQl" 
7.6 X 10^ 



2.4 X 10-1 
1.1 X 10^1 



1.1 X 10-" 
1.4 X 10-" 



stellar density profiles. The Sersic index of Aq-C's bulge is 
1.2 and that of Aq-D's bulge is 1.4. These values suggest 
that both bulges are 'pseudobulges'. The edge-on surface 
stellar density maps confirni that these bulges indeed 
show the characteristics of pseudobulges, they particularly 
resemble 'disky pseudobulges', while Aq-C also has features 
of a 'boxy bulge'. The b oxiness is consiste nt with the 
existence of a bar in Aq-C (|Athanassoulall2005l ) . 



~ 0.5 at a maximum; iFisher et al.l (|2009l ) reported that 
such galaxies have B/T < 0.6 in near-IR; All high-mass 
(M» ^ IO^^Mq) spirals whose bul ge indices are smalle r 
than 2 have _H"-band B/T < 0.4 (jweinzirl et all [200i). 
In our simulations, the values of H-hand BIT are 0.46 
and 0.21 for Aq-C and Aq-D, respectiveljy. Therefore 
Aq-C's B/T value is on the very high end of what has 
been observed and it might be out of range because we 
ignore the dust extinction which reduces the disc luminosity. 

Next, we compare the specific star formation rate in th e 
bulges with observed values. Following iFisher et al.l (|2009|), 
the bulge radius, J?xs, is defined as the radius at which the 
surface stellar density profile of a bulge is 25 % higher than 
that of a disc. The bulge mass, Mxs, is then given as: 



Mxs = M,tar(< i?xs) - 27rEa 



flxs 



-t/R, 



dr. (2) 



We list the bulge mass, the bulge radius, the star formation 
rate internal to the bulge, Mxs, and the specific star 
formation rate of the bulge, Mxs /Mxs, in Table [1] The 
specific star formation rates, 1.1 x lO^^^ yr^^ (Aq-C) and 
1.4 x 10"^"*^ yr~^ (Aq-D), are in good agreement wi th those 
for the inactive pseudobulges (JFisher et al.l [2009|). Pseu- 
dobulge formation via high-redshift starbursts provides a 
good explanation for the formation of inactive pseudobulges. 

Bulges of nearby galaxies ar e however less massive 
than 10^" Mq (|Fisher et al.ll2009l '). Aq-C's bulge (Mxs = 
2.2 X 10^°) seems to be too massive. This is because that 
the simulated galaxies have too much star formation at high 
redshift. Aq-C has formed half of th e final s t ellar mass by 
redshift 3.3 and Aq-D by redshift l.S. lLeitneJ (|2012l ') argued 
that current star-forming galaxies formed most of their stel- 
lar mass at redshift lower than 2. Results from abundance 
matching also show that Aq-C forms to o many stars at red- 
shift higher than 2 (|Moster et al.ll2012l . their Fig. 12) for its 
halo mass. The disagreement at high redshift indicates pos- 
sible failure of our galaxy formation model. We will discuss 
effects of this early star formation in Section [S] 



5 SUMMARY AND DISCUSSION 

We have performed A^-body/hydrodynamic simulations of 
galaxy formation in two Milky Way-sized haloes in a ACDM 
universe. The haloes are taken from the Aquarius project 
(|Springel et al.ll2003 ): Aq-C and Aq-D. Both galaxies form 
well defined discs; one has the bulge-to-total mass ratio, 
B/T ~ 0.6, and the other, B/T ~ 0.3, based on the surface 



* B/T becomes smaller in bluer bands (see also Table. lAH 



The formation histories of these pseudobulges are 
rather different from the standard picture, in which pseu- 
dobulges are believed to form through secular evolution of 
galaxy discs. In our simulations, the pseudobulges mainly 
form by high-redshift starbursts before redshift 2. The 
evolution of the surface stellar density profiles reveals that 
the pseudobulge components are already in place at redshift 
2-3 as disky components with small scale length. The mass 
of these central disky components at redshift 2 accounts for 
~ 70% and ~ 87% of the final pseudobulge mass of Aq-C 
and Aq-D, respectively. These progenitors of t he pseudob- 
ulges would be observed as high-redshift discs (jStark et al.l 
l200a ). In fact, when the galaxies have peak star formation 
rates, the star formation occurs in unperturbed gas discs. 
The high-redshift discs continuously change the orientation 
by exchanging the angular momentum with newly-accreted 
materials whose angular momentum is misaligned with the 
preexisting discs. In the course of the frequent flips, the 
discs are heated and turn into pseudobulge-like objects. 

The main stage of the formation of the pseudobulges 
ends soon after the starbursts by gas consumption due 
to the starbursts and the energetic winds that removes 
star-forming gas. When the host haloes are well established, 
the gas cooled in the haloes has larger angular momentum 
and gradually develops the discs with larger scale length 
than the pseudobulges. The secular evolution of the main 
discs that is mainly driven by bar instability in Aq-C and 
clumpy star formation in Aq-D further increases the mass of 
the pseudobulges. The secular evolution however accounts 
for only ~ 30% and ~ 13% of the mass of the pseudobulges 
of Aq-C and Aq-D, respectively. We therefore conclude 
that the main channel of the pseudobulge formation in our 
simulations is not the secular evolution of the discs. 

The specific star formation rates of the simulated 
bulges are co nsistent with thos e of the observed inactive 
pseudobulge (jFisher et al.l 12003 ). implying that such pseu- 
dobulges have formed via high-redshift starbursts. The 
bulge mass, Aq-C's bulge in particular, is however too large 
to compare with the observed counterparts. Aq-C indeed 
shows too much e arly star formation before redshift 2 
(jMoster et al.ll2012l '). We are thus likely to overestimate the 
importance of the high-redshift starbursts for pseudobulge 
formation especially for Aq-C. 

There seems to be two reasons why the simulated 
galaxies do not form classical bulges. One is that the host 
hal oes have rel a tively quiet merger histories as described 
in IWang et al.l (|201ll ). Consequently, the galaxies have 
not undergone mergers with mass ratios greater than 0.1 
since redshift 4. The merging mass ratios of the galaxies 
becomes smaller than those of their host haloes owing to 
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the feedback which is more effective in smaller galaxies. 
The other reason is that when the galaxies are quite gas 
rich when they undergo major mergers at high redshift. 
This is again due to our feedback that suppresses star 
formation very efficiently when the main progenitors are 
small. Wet (gas rich) mer gers are not as des tructive as 
dry (collisionless) mergers (IHopkins et al. 2010l) and often 
form discs as merger remnants (jSpringel fc Hernguistl l2005l : 
[Robertson et al.1 l2006l : iGovernato et all 120091 '). Anv suc- 
cessful simulations of d isc galaxy fo r mation must suppress 
early star formation IIBrook et al.l 2004 : lOkamoto et al.l 



l2005l : iGuedes et al.1 [20111 : IScannapieco et al.ll2012l '): we thus 
speculate that gas richness in high-redshift progenitors is 
the key to explain the existence of inactive pseudobulges in 
large disc galaxies. 

The formation scenario of pseudobulges by high-redshift 
starbursts presented in this paper provides a good expla- 
nation for the existence of inactive pseudobulges in early- 
type disc galaxies such as SO and S a galaxies. Pseudobulges 
do ex ist in early-type disc galaxies (jKormendv fc KennicuttI 
[2OOJ) and secular evolution may take too long to form such 
large pseudobulges. In our simulations, large inactive pseu- 
dobulges with the bulge-to-total mass ratios, B/T ~ 0.6 
and ~ 0.3; --^ 70 % and ~ 87 % of them have formed 
through high-redshift starbursts, respectively, and secular 
evolution adds the rest. A simulation of a smaller galaxies 
-^gai ~ 10^° Mq, for which feedback is more efficient than 
in Milky Way-sized galaxies, show ed that its pseudo bulge 
has formed by the secular evolution iBrook et al.l (|2012l ). We 
hence expect that different formation processes of pseudob- 
ulge operate depending on galaxy mass. 
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APPENDIX A: PROPERTIES OF THE 
SIMULATED GALAXIES AT REDSHIFT 

We here give some of the key properties of our simulated 
galaxies at redshift 0. In Table lAll we show the mass and 
i-band luminosity of them. We also present the best fit 
values of the parameters in the surface stellar density profile 
fitting. For reference, we also show the i-band bulge-to-total 
luminosity ratios, {B/T)i, which take lower values than the 
mass ratios because the discs contains more young stars 
than the bulges. The disc scale lengths seem to be too large 
compared with a typical value with observed galaxies of 
similar types (a few kpc). Such morphological details are 
quite sensitive to the implementation of star formation and 
feedback and beyond the scope of this paper. 
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Figure Al. Circular velocity profiles at redshift 0. The circular 
velocity profiles, Vc{r) = ^GM{< r)/r, are shown by the black 
solid lines. The blue dotted, red dashed, cyan dot-dashed, and 
yellow dot-dot-dot-dashod linos respectively indicate the contri- 
bution from the dark matter, stellar, gas, and baryonic (stellar 
plus gas) components. 
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Figure A2. Star formation efficiency as a function of halo mass. 
The vertical axis shows how many baryons in the halo locked 
into stars, where A/b = (f2b/^o)-'^200 a-rid M200 is the mass of a 
sphere whose mean densit y is 200 times the critical density. The 
black curve is taken from lGuo et al.l 1 I2OIOI ). which indicates the 
values required if a ACDM universe is to fit the observed SDSS 
stellar mass function. The red plus sign and blue asterisk indicate 
the values for Aq-C and Aq-D, respectively. 



In Fig. lAll we show the circular velocity profiles. 
Note that the circular velocity, Vc = \/GM{< r)/r, is 
different from the rotation speed presented in Fig. (31 We 
find that the circular velocity profiles are very flat and 
the values are consistent with the Milky Way's value 
(^ 200 km s~^;ISofue et al.ll2009l . see also G3-TO model in 



IScannapieco et al.ll2012f ). The stellar component dominates 



within ~ 4 kpc. Therefore within the bulges, stars dominate 
in mass and the outer region, where the discs exist, is 
dominated by dark matter. 



In Fig. IA2I we sho w how many ba r yons i n the haloes 
are lo cked into stars. Behroozi et al.l (|2010l ). IGuo et al.l 
(J2010h and lMoster et al.1 (|201Gh have pointed ou t only 20% 
of bar yons in a halo are locked into stars at most. IGuo et al.l 
(|2010|) showed the great majorit y of simulations lock to o 
many baryons into stars (see also I Scannapieco et al.ll2012l ). 
The values for out simulated galaxies are slightly low e r than 
that obtained for the SDSS galaxies by I Guo et al.l (|201G| ) 
but they are broadly consistent with their result consid- 
ering the fact that there must be significant halo-to-halo 
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Table Al. Mass and morphology of the simulated galaxies at redshift 0. The stellar mass, Mstar, and the i-band luminosity, Mi, of the 
galaxies are computed by using the stars within the galaxy radius. We present the best fit values of the parameters for the Scrsic bulges 
and explanation discs: So, Re and n in equation ((TJ, the disc central surface stellar density, T,^, and the disc scale length, R^. We also 
show the i-band bulge-to-total luminosity ratio, {B/T)i. 



Galaxy 


Mstar (Mq) 


Mi (mag) 


Ee (M0 kpc-2) 


Re (kpc) 


n 


So (Mq kpc-2) 


Rd (kpc) 


(B/T), 


Aq-C 
Aq-D 


4.0 X IQio 

3.1 X IQio 


-21.3 

-21.2 


2.2 X 10^ 
7.2 X 10* 


0.98 
1.1 


1.2 
1.4 


1.2 X 10* 
1.1 X 10* 
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0.38 
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Figure Bl. Star formation and merging histories of the simu- 
lated galaxies. The left and right panels show evolution of Aq-C 
and Aq-D, respectively. Upper panels: Star formation histories of 
the stars within the galaxy radius (O.lrvir) at redshift are shown 
by the black lines. The red vertical lines indicate merging mass 
ratios. All the mergers whose mass ratios are greater than 0.01 
are presented. Note that the number of outputs of Aq-D simu- 
lation is as four times large as that of Aq-C simulation; hence 
the number of the vertical lines is larger for Aq-D. Lower panels: 
Mass evolution of the main progenitor galaxies. The total (stellar 
-|- ISM), star-forming gas, and ISM masses in the main progeni- 
tor of each galaxy are shown by the black solid, blue dashed, and 
cyan dotted lines, respectively. 
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Figure B2. Mass evolution of the main progenitor haloes. Up- 
per panels: The mass content in the main progenitor haloes as a 
function of redshift. The black solid, yellow dotted, blue dashed, 
cyan dot-dashed, and red dot-dot-dot-dashed lines indicate mass 
in dark matter, stars, star-forming gas, warm gas, and hot gas, 
respectively. Lower panels: Evolution of the baryon fractions in 
the main progenitor haloes. The baryon fractions are normalised 
by the cosmic mean baryon fraction, Qi^/Qq. 



variation. The galaxy formation efficiency in Aq-D is lower 
than that in Aq-C since the formation redshift of Aq-D's 
host halo is more recent and the depth of its potential well 
is shallower; therefore the feedback is more efficient for 
Aq-D. 



APPENDIX B: ROLES OF INFLOWS AND 
OUTFLOWS 

Here we investigate the roles of inflows and outflows. First, 
we show the star formation and merging histories of the 
simulated galaxies in Fig. IBll We show the star formation 
histories of stars within the galaxy radius at redshift and 
the merging mass ratios. We also present the mass evolution 
of the galaxies and the mass in the ISM and star-forming 
gas in the lower panels. The mass evolution indicates that 



the galaxies are quite gas rich until the first starbursts at 
redshift ~ 5 for Aq-C and at redshift ~ 3.5 for Aq-D since 
our feedback effectively suppresses star formation in small 
galaxies. The major mergers at high redshift are thus wet. 

Roughly speaking, the starbursts above redshift 2 estab- 
lish the pseudobulges (see Fig.|6]) and therefore the quiescent 
star formation below redshift 2 seen in the star formation 
histories of the galaxies as a whole (Fig. IB1|I builds up the 
discs. Since mergers are not correlate with the star forma- 
tion activity, starbursts seem to be triggered by the rapid 
supply of low angular momentum gas. Indeed, the galaxy 
masses are dominated by the ISM when they grow rapidly 
as shown in the lower panels of Fig. IBll This suggests that 
gas accretion is much faster than star formation. 

Next, we illustrate the mass evolution of the main 
progenitor haloes of Aq-C and Aq-D in Fig. IB2I In the 
upper panels, we show the mass in dark matter, stars, 
star-forming gas, warm gas, and hot gas as functions of 
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redshift, where the star-forming gas is high density gas 
(tih > Tith = 1.6 cc~^), hot gas is that with T > 10^ K or 
that in the winds, and we define the rest as warm gas. 

Both haloes have virial mass of ~ 10^^ Mq at redshift 
as we chose the haloes of the Milky Way's halo mass. 
The amount of the star-forming gas becomes largest at 
redshift ~ 5 in Aq-C and at redshift ^ 3.5 in Aq-D. 
Below these redshifts the mass of the star-forming gas 
monotonically decreases with decreasing redshift owing to 
the star formation and feedback. 

The effect of winds is depicted in the lower panels of 
Fig. IB2I in which we plot the baryon fraction in the main 
progenitor haloes. The baryon fraction is normalised by 
the cosmic mean value il^/ilo, that is, if a halo has not 
lost any baryons, the value becomes unity. The baryon 
fraction is high at high redshift and decreases towards lower 
redshift owing to the winds. There is an interesting feature 
between redshift 5-9 in Aq-C and 3.5-7 in Aq-D, where 
the baryon fractions increase towards lower redshift. At 
these redshifts, the haloes glow rapidly and the gas inflow 
competes the outflow. Consequently, the baryon fractions 
are peaked at redshift ~ 5 in Aq-C and at redshifts ~ 3.5 
in Aq-D. These redshifts coincide with the redshifts at 
which the mass of the star- forming gas becomes largest; the 
starbursts occur at these epochs. The baryon fractions then 
shaxply drop because strong supernova feedback following 
the starbursts blow out the gas from the haloes. During 
the quiescent star formation (redshift < 2), the baryon 
fractions keep decreasing and they are only one third of 
the cosmic mean at redshift 0. After the first starbursts, 
the total baryon mass in the haloes are almost constant, 
while the baryon fractions decreases. This means that the 
mass flux in outflows balances that in the inflows. The 
outflows preferen t ially r emove low angular momentum gas 
l|Governato et all l2010l : iBrook et ahl 120111 ): on the other 
hand, low-redshift inflows in general bring in high angular 
momentum gas. 



Aq-C 



Aq-D 



Fig. IB3I shows the distribution of angular momentum 
of inflows and outflows during the starbursts. Outflows (in- 
flows) are defined as gas moving radially outward (inward) 
between 0.3rvir and 0.7rvir where wind particles have already 
been recoupled to the hydrodynamic interactions. The z- 
direction is chosen as the direction of the spin of the galaxy 
at a given redshift. We find that the angular momentum dis- 
tributions of the outfiows are always peaked at jz — show- 
ing that the outfiows consist of low angular momentum gas. 
At redshift < 3, the distribution of the angular momentum 
of the infiows show peaks at j^ > 0. The outfiows however 
remove a large amount of high angular momentum gas too. 
Feedback and star formation prescriptions that remove low 
angular mo mentum gas more efficiently would yiel d smaller 
bulges (see lGuedes et al.ll201ll : lMaccio et al.ll2012l ). 



APPENDIX C: NUMERICAL CONVERGENCE 

In order to conduct convergence studies, we also simulated 
the Aq-C halo at two lower resolutions with particle masses 
~ 8 and ~ 64 times larger. Here we employ the same 



a 0.0003 



0.0000 I i""T 



0.0002 
0.0001 : 
0.0000 



outflow 
inflow 




z = 2.9 







-10000-5000 5000 10000 
j, [kpc km s-'] 



-10000-5000 5000 10000 
j, [kpc km s-'] 



Figure B3. Distribution of angular momentum of outflows and 
inflows. Outflows (inflows) are defined as gas moving radially out- 
ward (inward) between O.Srvir and 0.7rvir. The red solid and blue 
dotted lines indicate outflows and inflows, respectively. From top 
to bottom, the distributions at redshift ~ 4, 3, and 2 are shown. 



Table CI. Numerical parameters employed for the three differ- 
ent resolution simulations: dark matter and SPH particle masses 
and the gravitational softening length in physical units. The grav- 
itational softening length is kept in constant in comoving units 
until redshift 3 and then fixed in physical units, i.e. the values 
shown in the table. 





A/dm (Mq) 


MsPH (Mq) 


ephys (pc) 


Aq-C-4 
Aq-C-5 
Aq-C-6 


2.6 X 10^ 
2.1 X 106 

1.7 X 10^ 


5.8 X 10* 
4.7 X 10^ 
1.7 X 10^ 


257 
514 
1028 



naming convention as Aquarius, labelling the three runs (in 
order of decreasing resolution) Aq-C-4, Aq-C-5 and Aq-C-6. 
Aq-C-4 is the one we have called 'Aq-C in the main text of 
the paper. Table [CTI lists the numerical parameters of each 
simulation. 

The only physical parameter that depends on the 
numerical resolution is the star formation threshold density, 
nth. The star formation threshold density for Aq-C-4, -5, 
and -6 are nth = 0.1, 0.4, and 1.6 cc~^, respectively. A 
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Figure CI. Formation histories of stars witliin galaxy radius in 
three different resolution simulations. The upper panel shows the 
differential star formation histories while the lower panel presents 
them in cumulative form. The black solid, blue dotted, and red 
dashed lines corresponds to Aq-C-4, Aq-C-5, and Aq-C-6, respec- 
tively. 



factor of 4 higher density is adopted for a factor of 8 higher 
mass resolution. This scaling is chosen since, for irradiated 
nearly isothermal gas with an isothermal density profile, 
halving the gravitational softening length will increase the 
maximum density that is resolved by a facto of 4. 

Using the same set of the simulations, convergence 
studies on t he sa tellite population was carried out by 
IParrv et al.l (|2012f ) and they found reasonably good 
convergence of various key properties of the simulated 
satellite galaxies; such as star formation histories, mass 
of the satellites and luminosity functions. Comprehensive 
co nvergence studies on the central galaxy were presented 
in IScannapieco et al.l ()2012l ') but only Aq-C-5 and Aq-C-6 
were used. 

In Fig. ICll we show the star formation histories of 
these galaxies. The results from three different simulations 
agree quite well, although there are some differences in the 
differential star formation histories during the starbursts. 
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Figure C2. Mass- weighted normalised distribution of orbital 
circularity, jz/jc{E), of all stars lying within Tvi,- of the centre, 
excluding stars identified as members of satellite galaxies. Aq-C- 
4, Aq-C-5, and Aq-C-6 are indicated by black solid, blue dotted, 
and red dashed lines, respectively. For reference, we also show 
the result for Aq-D (corresponding to the level 4 resolution) by 
the green dot-dashed line. A peak at jz/jc(E) ~ 1 indicates the 
existence of a disc component. 



Such level of difference will occur even if we only change a 
seed for a random number sequence because of stochasticity 
of the star formation. 

In order to compare the morphologies of the galaxies, 
we analyse the distribution of the orbital circularity, defined 
below, of stars in the host haloes. We exclude stars be- 
longing to the satellites. For each star lying within rvir, we 
compute jz, the component of specific angular momentum 
parallel to the 2-axis. We then compute the specific angular 
momentum, jc{E), of a prograde circular orbit with the 
same binding energy as the particle in question. The ratio 
.iz/.ic(E) defines the orbital circularity (jAbadi et al.l 120031 : 



lOkamoto et al.ll200a). The way to compute this quantity is 
described in lOkamoto et al.l (|2010l ) in detail. 



Fig:. IC2l shows the normalised distributions of the orbital 
circularity of stars within rvir for three difference resolution 
simulations. We also present that for Aq-D for reference. A 
cold disc component has jz/jc{E) ~ 1 and such a compo- 
nent is evident in all galaxies. There is however no obvious 
convergence feature. For example, the fraction of stars with 
retrograde orbits (i.e. j^ < 0) is largest in Aq-C-5; there- 
fore Aq-C-5 is the most spheroid-dominated galaxy of the 
threqj. The position of the second highest peak in Aq-C- 
4 ijz/jc — 0.5) is different from that in lower resolution 

^ If we a ssume the n e t rota tion of a spheroidal component is 
zero as in lAbadi et al.l 112003) , the mass of a spheroid is given by 
2MUz < 0). 
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counterparts, Aq-C-5 and Aq-C-6 {jz/jc = 0). Despite the 
fact that there seems to be no convergence feature, it is en- 
couraging that changing resolution does not totally change 
the morphology, that is, they all form well-defined discs. By 
comparing Aq-C-4 with Aq-D, we find Aq-D is more disc 
dominated than Aq-C-4. This is consistent with the result 
by the decomposition based on the surface stellar density 
profiles. 

This paper has been typeset from a TJjX/ Wlff^ file prepared 
by the author. 



